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Abstract 

The method of self-similar factor approximants is completed by defining the ap- 
proximants of odd orders, constructed from the power series with the largest term of 
an odd power. It is shown that the method provides good approximations for tran- 
scendental functions. In some cases, just a few terms in a power series make it possible 
to reconstruct a transcendental function exactly. Numerical convergence of the factor 
approximants is checked for several examples. A special attention is paid to the pos- 
sibility of extrapolating the behavior of functions, with arguments tending to infinity, 
from the related asymptotic series at small arguments. Applications of the method 
are thoroughly illustrated by the examples of several functions, nonlinear differential 
equations, and anharmonic models. 



PACS numbers: 02.30.Lt; 02.30.Mv; 03.65.Ge 
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1 Introduction 



The problem of reconstructing functions from their asymptotic series is widely known to be 
of great importance for a variety of applications [1]. Probably the most often used technique 
allowing for such a reconstruction is that based on Pade approximants [2], though there also 
exist other more complicated methods [3]. 

Recently, a novel approach has been advanced for reconstructing functions from the 
related asymptotic series, called the method of self-similar factor approximants [4-6]. The 
form of these approximants is derived by means of the self-similar approximation theory 
[7-15]. The factor approximants were shown [4-6] to be more general and more accurate 
than Pade approximants. However, several points in the method of factor approximants 
[4-6] have not been investigated. 

The most important point, which has remained unclear, is how to construct the factor 
approximants of odd orders'! The problem is that the standard procedure [4-6] requires that 
the generic series, the factor approximants are built from, be of even order. So that the 
large amount of information, contained in the odd orders of the series, could not be properly 
used. Here we advance a uniform approach for treating both odd as well as even orders 
of power series and we demonstrate, by a number of examples, good performance of this 
general method. 

Another, rather nontrivial, question is why the factor approximants could provide high 
accuracy for transcendental functions, obtained from their asymptotic series. We give an 
explanation for this effect, which is unique for the resummation methods based on asymptotic 
series. Recall that Pade approximants, having the structure of rational functions, are usually 
not so good for approximating transcendental functions. Moreover, we show that some 
transcendental functions can be reconstructed, by means of the factor approximants, exactly 
from a few terms of an asymptotic series. 

In the previous papers [4-6] on the method of self-similar factor approximants, we mainly 
considered the low-order approximants of even order, because of which their numerical con- 
vergence could not be properly analyzed. Now, we possess both even and odd orders of 
these approximants, and we study their numerical convergence for various cases, calculating 
high-order approximants, up to 20-th order. 

We pay a special attention to the possibility of extrapolating physically motivated func- 
tions to the whole range of their variables. We demonstrate that the method of self-similar 
factor approximants can serve as a tool for defining the values of functions at infinity from 
their asymptotic expressions near zero. 

Summarizing, the main results of the present paper are as follows: 

(i) The method of self-similar factor approximants is extended by constructing the factor 
approximants of odd orders, in addition to those of even orders. 

(ii) It is shown that transcendental functions can be well approximated. A unique feature 
of the method is that some transcendental functions can be reconstructed exactly. 

(iii) The method allows for the extrapolation of solutions to nonlinear differential equa- 
tions from asymptotically small to finite variables. 

(iv) The characteristics of quantum anharmonic models can be extrapolated from the 
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region of asymptotically small coupling parameter to the whole range of the latter, including 
the case of the coupling parameter tending to infinity. 

(v) For all considered examples, the factor approximants of high orders are calculated, 
demonstrating the existence of numerical convergence. 



2 Self-similar factor approximants 

Let us be interested in defining a real function f(x) of a real variable x G M. The extension 
of the method to complex- valued functions of complex variables is also possible, but, first, 
let us consider a slightly simpler case of real functions and variables. Suppose, we know only 
the behavior of the function at asymptotically small values of x — > 0, where we can get the 
sequence {fk( x )} of the expansions 

k 

fk(x) = a nX H , (!) 
n=0 

with k — 0,1,2, . . .. The sequence, generally, can be divergent. It is known [1] that a function 
f(x), analytic in the vicinity of x — 0, uniquely defines its asymptotic series. The converse is 
not always true. But in what follows, we consider the situation, when there is a one-to-one 
correspondence between the function and its asymptotic series, since solely then the problem 
of reconstructing functions from their expansions acquires sense. 

Without the loss of generality, we may assume that in expansion (1), 

«o = /*(0) = f (x) = 1 . (2) 
Really, if instead of expansion (1), we would have 

f k \x)=f^(x)J2a n x n , 

n=0 

with a given function f(°\x), then we could immediately return to Eq. (1) defining 

f (k \x) 



A 0*0 



/ (0) (^)«o ' 



The self-similar factor approximants of even orders k — 2p — 2, 4, 6, . . . are given [4-6] 
by the form 

= n (i + > ( 3 ) 

i=i 

whose parameters Ai and n» are obtained from the re-expansion procedure, that is, by 
expanding Eq. (3) in powers of x and comparing the results with the given expansion 
(1), equating the like-order terms. This accuracy-through-order procedure yields the set of 
equations 

Y,mAi = B n (n = l,2,...,2p), (4) 

i=i 
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with the right-hand sides 



(_-i\n-l Jn 

B n = y^— Km — \nf k (x). (5) 
(n — 1)! x->o dx n 

In each approximation, the parameters Ai, rii, and B n , of course, depend on the approxima- 
tion number k = 2p. However, in order to avoid too cumbersome notation, we do not mark 
explicitly this dependence, which is assumed to be evident. 

Each factor in form (3) contains two parameters, Ai and rii. This is why these approx- 
imants could be straightforwardly defined only for the even-order expansions f2 P (x). But 
how should we proceed having odd-order expansions (1) with k — 2p + 1 = 1, 3, 5, . . .? For 
the latter, we could write the form 

p+1 

/2 P +iW = n( i +^r , (6) 

i=i 

with the parameters satisfying the set of equations 

p+1 

Y,rnA? = B n (n = l,2,...,2p+l) . (7) 

i=i 

But the problem is that form (6) contains 2p + 2 unknown Ai and rii, while only 2p + 1 
equations of set (7) are available. One equation is lacking. 

To overcome this problem, let us notice that expression (6) is invariant under the scaling 
transformation 

x — > Xx , Ai -> X^Ai . (8) 

Then, taking A — > A^ 1 and using the renotation Ai/Ai — > A h we come to the same form 
(6), but with 

Ai = 1 (fc = 2p+l) . (9) 

Complementing the set of 2p + 1 equations (7) by the scaling condition (9), we get 2p + 2 
equations for 2p + 2 unknowns. In this way, we now can construct the factor approximants 
of odd orders. 

If the sought function f(x) can be reduced to the factor form f%(x), then it is obvious 
that this function can be reconstructed exactly for all approximants, starting with the given 
order k, since both f(x) and f^{x) yield the same series of the fc-th order. Thus, there exists 
a class of functions , having the form of f^(x), which can be reconstructed exactly 

from their asymptotic series by means of the factor approximants. This class of exactly 
reproduceable functions can be noticeably extended by including those functions that can 
be defined as limits of the form fZ(x) with respect to the values of some Ai and n» tending to 
either zero or infinity. Keeping in mind such limits, we shall denote them as lim^.^.}. All 
functions, which can either be reduced to the factor form fZ(x) or obtained from the latter 
by means of limiting procedures with respect to their parameters, form the factor class 

F=\f;(x) lirn r k (x)\ . (10) 

{Ai,rii} I 



4 



As is clear, a function can be exactly reproduced by the factor approximants if and only if 
it pertains to the factor class (10). This is a rather wide class, including rational as well as 
irrational functions. In the following section, we show that some transcendental functions 
can also be reproduced exactly. 

3 Reconstruction of transcendental functions 

Consider such an entire transcendental function as the exponential 

oo T n 

m = = e -7 • (ii) 



n=0 Ul 

From the expansion f 2 (x) — 1 + x + \x 2 , we find the factor approximant 

f*(x) = lim(l + ax) 1/a = e x . (12) 

Similarly, from the form fs(x) at x — > 0, we obtain 

f£(x) = lim(l + x) b/(1 ' b) (l + bx) 1/h{l - b) = e x . (13) 

This procedure can be continued to any order, showing that the transcendental function e x 
is exactly reproduceable by the factor approximants of any finite order k > 2, 

f* k (x)=e* (k>2). (14) 

The hyperbolic functions, such as sinhx and coshrr, can also be reconstructed exactly, if 
one notices that the series, corresponding to 2sinha; = e x — e~ x and 2coshx = e x + e~ x , can 
be represented as the sums of two series for the exponentials, which, as is shown above, are 
reproduceable exactly. The trigonometric functions sinx and cos a; can also be reproduced 
exactly, if one treats them as the analytic continuations of the related hyperbolic functions 
to the imaginary axis, when x is replaced by ix. 

The fact that some transcendental functions can be reproduced by factor approximants 
exactly does not imply that all transcendental functions are exactly reproduceable. How- 
ever, the very existence of exactly reproduceable transcendental functions explains why other 
transcendental functions, though not exactly reproduceable, nevertheless, can be well ap- 
proximated by the factor approximants. 

As an example, let us consider the series for sin a;, 

oo 



(-1)" 

n =o (2n-l)! 



sinx = J2 7KZ — t x2n+1 ■ (I 5 ) 



Note that for the alternating series, it is convenient to shift them by a constant, so that all 
factor approximants for the shifted function be sign defined, removing this shift afterwards. 
The actual value of the shifting constant plays no role, provided the considered function 
becomes sign defined. In the present case, we, first, consider sinx + const, and remove the 
constant at the end. It is important to emphasize that the factor approximants, obtained 
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from the Taylor series, are much more accurate than the latter. For an illustration, we show 
in Fig. 1 the factor approximant of the order k = 18, derived from series (15), and the 
Taylor series of the same order k = 18, as compared to the exact function sin a;. As is seen, 
the factor approximant is essentially more accurate. 

Constructing the factor approximants for the tangent, 

oo o2n i 

tan* = £(-l) B 2 2 « — — - B 2n x 2n ~ l , (16) 

71=0 \ An )- 

where B n is the Bernoulli number, we get such a good accuracy for the approximants of 
the order k > 4, that they practically coincide with the exact function. Therefore, we 
concentrate on defining the point of singularity, which for tan* occurs at * = ±7r/2. For 
each factor approximant, we find the point of its singularity x k and compare it with the 
exact value n/2. We consider the positive value * > 0, since for * < 0, the behavior is 
similar, being antisymmetric. The accuracy of the approximants can be characterized by the 
percentage error 

e k = • 100% . (17) 

7T/ Z 

The error quickly diminishes with the approximation order, demonstrating fast numerical 
convergence, for instance, e 2 = —6.8%, e 3 = —5.4%, e A = 0.13%, e 5 = 0.096%, and e 6 = 
-0.00035%. 

Numerical convergence is well seen for different transcendental functions. We have con- 
sidered a variety of the latter. For illustration, we show in Fig. 2 how the factor approximants 
approach the function 

/(*) = x + cos* . (18) 

Often transcendental functions appear as solutions of differential equations. Let the 
function y(t) of time t > be defined by a differential equation. The simplest way to 
construct its solutions in the form of the factor approximants is as follows. We may look for 
the solution at small t — > 0, representing it as a series 

n 

And then, it is straightforward to derive the factor approximants from the given series, as is 
explained in Sec. 2. 

Let us consider the nonlinear ordinary differential equation 

% + V + V> = 0, (19) 
describing a nonlinear oscillator with the cubic potential 

The solution to this equation (19) corresponds to the oscillations around the center (y = 
0, y — 0), where the overdot, as usual, means time derivative, during the time interval 
< t < t . At the point t Q , depending on initial conditions, the solution diverges as 
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y oc (t — to) -2 - Constructing the solutions to this equation from the related time series 
at t — > 0, as is explained above, we obtain the factor approximants y%(t) providing the 
approximate solutions for longer times. In Fig. 3, the factor approximants yl 4 {t) and yl e {t) 
are shown in comparison with the exact numerical solution of Eq. (19) with the initial 
conditions y(0) = and y(0) = I. For these initial conditions, the divergence occurs at 
t = 6.5. We have analyzed the behavior of high-order approximants yl(t), with k > 15, 
finding that they provide a good extrapolation from the region of asymptotically small times 
t — > to the temporal interval < t < 5. 

As another example of nonlinear differential equations, let us take the Rayleigh equation 

d 2 V . dy e (dy\ 3 

describing an oscillator perturbed by the right-hand side of Eq. (20), with < £ -C 1. 
Written in the normal form, Eq. (20) represents a two-dimensional dynamical system having 
one unstable fixed point (y — 0, y — 0), with the Jacobian eigenvalues = (e± \J e 1 — 4)/2. 
By Poincare-Bendixon theorem [16], there exists a limit cycle around this unstable fixed 
point. We look for the solution of Eq. (20), with e = 0.1, under the initial conditions 
2/(0) = 0, y(0) = 1. As is described above, we, first, find the power-series solution at t — > 0, 
from which we construct the factor approximants yl(t). In Fig. 4, the approximants yl 8 (t) 
and ylg(t) are compared to the exact numerical solution of Eq. (20). Again we observe 
that the factor approximants of the order k > 19 well extrapolate the time-series solution at 
t — > to the finite region < t < 6. 

It is worth noting that the problem of extrapolating time series is of high importance 
for the possibility of predicting the evolution of different complex systems such as markets 
and societies [17]. The examples, considered above, demonstrate that the self-similar factor 
approximants provide us with the tool of extrapolating the time series for small times t — > 
to the region of finite times, that is, making it possible to give forecasts, at least, for some 
near future. 



4 Extrapolation for anharmonic models 

Many interacting physical systems have the structure analogous to anharmonic models. This 
is why the latter serve as a touchstone for testing different approximation theories. The 
question to be investigated in the present section is as follows. Suppose, we are able to find 
a function of interest only by means of perturbation theory with respect to the coupling 
parameter g. We keep in mind here the simple perturbation theory in powers of g, without 
invoking any optimization procedures (see review-type articles [12,15]). Then, as is well 
known, the sought function has the form of an asymptotic series 

/(j)^EM B (<7->0), (21) 

n 

which diverges for any finite g. Could we, being based solely on this asymptotic expansion, 
extrapolate by factor approximants the function to the whole region of the coupling param- 
eter g G [0, oo)? It is especially interesting whether we could predict the behavior of the 
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function at very large couplings g — > oo. The latter has, clearly, no sense for the power series 
derived for g — > 0. Pade approximants are also not of the help for this purpose, since the 
limit g — > oo for a Pade approximant P[m/n\{.q) is n °t defined, as 

P[M/N](g) oc g M ~ N (g->oo) 

can tend to infinity, zero, or a constant, depending on M and N, so that, as g — > oo, 

f oo, M > N 
P[M/N](g) -> < const, M = N 
{ 0, M < N . 

First, let us consider the so-called partition function of the zero-dimensional y? 4 model, 

Z(s) = 4= r e-^-^d V . (22) 

The power series (21), corresponding to Eq. (22), possess the coefficients 

(-l) n / 1 
^7™! V 2 

which factorially increase at large n, because of which series (21) diverges for any finite g. 

Being based on series (21), we construct, according to the general scheme of Sec. 2, the 
factor approximants Z%(g) and check their accuracy by calculating percentage errors 

e k (g) ss Zli9 l~ g) Zi9) ■ 100% (23) 

with respect to the exact value (22). The results for different approximation orders and 
coupling parameters are presented in Table 1. This shows that the even-order factor ap- 
proximants demonstrate the uniform numerical convergence for all g, and the errors of the 
odd-order approximants oscillate in the frame of the errors of the even-order approximants. 
At large g, the partition function (22) has the asymptotic behavior 

Z(g) ~ 1.022765 g~ 1/4 (g -> oo) . (24) 

It is interesting whether the factor approximants (g) , based solely on the power expansion 
(21) at small g — > 0, could, nevertheless, catch the behavior of Z(g) at large g — > oo. The 
most nontrivial would be to predict the noninteger index —1/4 in the asymptotic form (24). 
For the found factor approximants, we study the limiting expression 

Z* k (g) ~ c k g- a * (g - oo) (25) 

and define the percentage error of the predicted index a fe as 

. . oth — 0.25 „ 

= fc Q 25 • 100% . (26) 
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In Table 2, we present the parameters of the asymptotic form (25), characterizing the behav- 
ior of the factor approximants at large coupling g — > oo, together with the percentage error 
(26) for the index a k - Again, we see that the even-order approximants display numerical 
convergence, and the odd-order approximants oscillate around the values of the even-order 
ones. The index a k , as well as the amplitudes c k , predict the large-coupling limit (24) with 
an error of about 20%. This is not as bad, if we remember that all calculations are based 
solely on the weak-couping expansion (21), valid for g — > 0, and no additional optimization 
procedures have been involved. 

Another standard touchstone for testing new approximations is the calculation of the 
ground-state energy for the one-dimensional anharmonic oscillator with the Hamiltonian 

"—klZ + k*^' (27) 

in which x G (—00,00) and g G [0,oo). The weak-coupling expansion (21) for the ground- 
state energy E(g) is known [18,19]. In the latter references, one can find the corresponding 
values of the expansion coefficients a n . 

Following the standard procedure of Sec. 2, we construct the factor approximants E%(g) 
for the ground-state energy E(g) of the anharmonic oscillator with the Hamiltonian (27) and 
compare the obtained El(g) with the exact numerical data [20]. The percentage errors 

e k (g) = Eli9 l^ i9) ■ 100% (28) 

of the factor approximants are given in Table 3. Our results show that the even-order 

and odd-order approximants form two sequences {E£ p (g)} and {E% +1 (g)} each of which 
uniformly converges for all g. 

We again pay a special attention to the ability of factor approximants to predict the 
behavior of the sought function at infinitely large coupling parameter g — > 00. Then the 
ground-state energy behaves as 

E{g) ~ 0.667986 g 1/3 . (29) 

We analyze the factor approximants E%(g), obtained on the basis of the weak-coupling 
expansion (21) at g — > 0, in the strong-coupling limit g — > 00 and compare their asymptotic 
form 

E* k (g) ~ b k gh (g - 00) (30) 

with the exact asymptotic behavior (29). The results are presented in Table 4, where also 
the percentage error 

em = • 100% 

of the predicted index (3 k is shown. The even-order approximants demonstrate numerical 
convergence, and the errors of the odd-order approximants fluctuate around those of the 
even-order approximants. The errors of the amplitudes b k are close to 10% for all k > 2. 
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5 Conclusion 



The method of self-similar factor approximants has been completed by defining a general 
way for constructing both even-order as well as odd-order approximants. It is shown that 
transcendental functions can be well approximated by factor approximants. Moreover, some 
transcendental functions can be reconstructed exactly from just a few terms of their asymp- 
totic expansions at small arguments. For all cases considered, the even-order approximants 
demonstrate uniform numerical convergence. The accuracy of the odd-order approximants 
is close to that of the even-order ones. 

It is shown that the solutions to nonlinear differential equations can be extrapolated from 
the values of asymptotically small variables to the finite region of the latter. 

A special attention has been paid to the possibility of predicting the asymptotic behavior 
of functions at their arguments tending to infinity from the knowledge of their asymptotic 
expansions at zero. Such an extreme extrapolation is shown to be feasible by means of the 
factor approximants. 

Also, we have studied the problem of a possible evaluation of the accuracy of the factor 
approximants, when the exact solution would not be available. For this purpose, we have 
compared, for the cases analyzed above, the quantities (/£ — f£_ 1 )/2 and — /, where / 
is the known exact numerical solution. It turned out that these quantities are close to each 
other. Hence the difference (/£ — f^._ 1 )/2 can serve as a reasonable estimate for the accuracy 
of the fc-th order factor approximant /£. 
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Figure Captions 



Fig. 1. Factor approximant f* s (x) (dashed line) and the Taylor series of the same order 
k — 18 (dashed-dotted line), as compared to the exact sinx (solid line). 

Fig. 2. Factor approximants f${x) (dashed-dotted line), fg{x) (dotted line), and f* (x) 
(dashed line), as compared to the exact function (18) (solid line). 

Fig. 3. Solutions to the nonlinear oscillator equation (19), with the initial conditions 
y(0) = and y(0) = 1, obtained from the time series at t — > in the form of the factor 
approximants y%{t). Here shown are y*±(t) (dotted line), y* & (t) (dashed line), and the exact 
numerical solution (solid line). 

Fig. 4. Solutions to the Rayleigh equation (20) in the form of the factor approximants 
yl 8 (t) (dotted line) and y* 9 (t) (dashed line), compared to the exact numerical solution (solid 
line) . 



12 



Table Captions 



Table 1. Percentage errors of the factor approximants Z%(g) for the partition function 
of the anharmonic model (22). 

Table 2. Amplitudes c k and indices a k for the strong-coupling limit (25) of the factor 
approximants for the partition function (22), together with the percentage errors (26) for 
the predicted indices. 

Table 3. Percentage errors of the factor approximants E%(g) for the ground-state energy 
of the anharmonic oscillator. 

Table 4. Amplitudes b k and indices (3 k for the strong-coupling limit (30) of the factor 
approximants for the ground-state energy of the anharmonic oscillator, with the percentage 
errors e k (P) of the predicted indices (3 k . 
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Table 1 



k\g 


0.1 


1 


5 


10 


100 


2 


0.2 


5.4 


19 


27 


70 


3 


-0.07 


-4.6 


-17 


-23 


-41 


4 


0.02 


2.0 


9.4 


15 


42 


5 


-0.008 


-1.3 


-5.7 


-8.0 


-13 


6 


0.003 


0.9 


5.8 


9.7 


30 


7 


-0.001 


-0.5 


-2.0 


-2.5 


-1.6 


8 


0.0006 


0.5 


3.9 


6.9 


23 


9 


-0.0002 


-0.2 


-0.5 


-0.2 


3.8 


10 


0.0001 


0.3 


2.8 


5.2 


19 


11 


-0.00004 


-0.05 


0.2 


0.8 


6.4 


12 


0.00004 


0.2 


2.1 


4.1 


16 
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Table 2 



k 


Cfc 




e k (a) 


2 


0.823 


0.090 


-64 


3 


0.917 


0.346 


38 


4 


0.806 


0.129 


-48 


5 


0.878 


0.255 


2.0 


6 


0.806 


0.148 


-41 


7 


0.860 


0.223 


-11 


8 


0.810 


0.161 


-36 


9 


0.850 


0.209 


-16 


10 


0.814 


0.170 


-32 


11 


0.845 


0.202 


-19 


12 


0.819 


0.178 


-29 


13 


0.842 


0.199 


-20 


14 


0.824 


0.182 


-27 


15 


0.841 


0.197 


-21 


16 


0.828 


0.187 


-25 


17 


0.840 


0.196 


-22 
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Table 3 



k\g 


0.01 


0.3 


1 


200 


2 


-0.07 


-2.0 


-7.4 


-53 


3 


-0.07 


1.4 


9.0 


256 


4 


-0.07 


-0.4 


-2.4 


-35 


5 


-0.07 


0.2 


2.2 


49 


6 


-0.07 


-0.1 


-1.0 


-25 


7 


-0.07 


0.06 


0.7 


14 


8 


-0.07 


-0.03 


-0.5 


-19 


9 


-0.07 


0.02 


0.3 


2.6 


10 


-0.07 


-0.01 


-0.3 


-15 



16 



Table 4 



k 


h 


Pk 


eM 


2 


0.729 


0.176 


-47 


3 


0.611 


0.590 


77 


4 


0.755 


0.231 


-30 


5 


0.669 


0.409 


22 


6 


0.756 


0.257 


-23 


7 


0.696 


0.351 


5.3 


8 


0.752 


0.272 


-18 


9 


0.710 


0.328 


-1.7 


10 


0.748 


0.282 


-16 


11 


0.718 


0.317 


-4.8 


12 


0.743 


0.289 


-13 


13 


0.721 


0.312 


-6.3 


14 


0.739 


0.294 


-12 


15 


0.723 


0.309 


-7.2 


16 


0.736 


0.298 


-11 


17 


0.725 


0.308 


-7.5 
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Figure 1: Factor approximant fi 8 (x) (dashed line) and the Taylor series of the same order 
k — 18 (dashed-dotted line), as compared to the exact sinx (solid line). 
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COS X + X 




Figure 2: Factor approximants f£(x) (dashed-dotted line), fg(x) (dotted line), and ff (x) 
(dashed line), as compared to the exact function (18) (solid line). 
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Figure 3: Solutions to the nonlinear oscillator equation (19), with the initial conditions 
y(0) = and y(0) = 1, obtained from the time series at t — > in the form of the factor 
approximants y%{t). Here shown are yl 4 (t) (dotted line), y* 6 (t) (dashed line), and the exact 
numerical solution (solid line). 
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Figure 4: Solutions to the Rayleigh equation (20) in the form of the factor approximants 
yl 8 (t) (dotted line) and y* g (t) (dashed line), compared to the exact numerical solution (solid 
line) . 
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